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Abstract 

We examine the sensitivity of recent LHC searches to signatures of supersymmetry with 
i?-parity violation (RPV). Motivated by naturalness of the Higgs potential, which would favor 
light third-generation squarks, and the stringent LHC bounds on spectra in which the gluino 
or first- and second- generation squarks are light, we focus on scenarios dominated by the 
pair production of light stops. We consider the various possible direct and cascade decays of 
the stop that involve the trilinear RPV operators. We find that in many cases, the existing 
searches exclude stops in the natural mass range and beyond. However, typically there is 
little or no sensitivity to cases dominated by UDD operators or LQD operators involving taus. 
We propose several ideas for searches which could address the existing gaps in experimental 
coverage of these signals. 
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1 Introduction 

The very successful first years of running at the Large Hadron Collider (LHC) have created an 
exciting atmosphere for particle physics. With 5 fb~^ of 7 TeV data analyzed by ATLAS and 
CMS, the LHC experiments have significantly tightened bounds on models of new physics. For 
decades, supersymmetric (SUSY) theories have been some of the leading candidates to resolve the 
few tensions with the Standard Model (SM). The minimal supersymmetric standard model (MSSM) 
has been viewed as one of the most well- motivated candidates for TeV scale physics, so the LHC's 
overwhelming validation of the SM and the apparent absence of MSSM signals thus far has proven 
quite startling. At the same time, the current null results constrain parameter space, providing the 
community with a vision of where to search for the MSSM, if it is present. 

In the MSSM, obtaining the electroweak symmetry breaking scale without fine tuning requires 
that the SUSY breaking soft mass parameters of the right-handed stop and the left-handed stop- 
sbottom doublet are not too large. This condition, naturalness, suggests that the two stops and at 
least one of the sbottoms should not be heavier than roughly 500 GeV (for a recent review, see [1]). 
However, results from the LHC indicate that, in most models, gluinos and most of the squarks 
must be beyond this natural range of masses. Depending on the details of the spectrum, the lower 
bounds on the masses of these particles vary between about 600 GeV and above 1 TeV for both 
i?-parity conserving (see, e.g., [2-5, 1]) and i?-parity violating (RPV) scenarios (see, e.g., [6, 7]). 
This motivates us to focus our attention on scenarios in which the third generation squarks are 
much lighter than all other colored superpartners. 



1 



Additionally, the discovery of a particle consistent with a Higgs boson at m/j ^ 126 GeV suggests 
that the two stops should be split by more than a few hundred GeV, however this does not need 
to be the case in extensions of the MSSM (see, e.g., [8]). Extensions of the MSSM are strongly 
motivated by the fact that a Higgs of this mass makes the MSSM relatively fine-tuned (unnatural) 
for any choice of the soft parameters. Depending on the nature of the extension, production and 
decay processes may be affected. However, since no very compelling extension has been suggested 
thus far and it is conceivable that the processes we consider will not be affected, we believe that 
studying MSSM processes would be most useful at this stage. Without knowledge of whether the 
heaviness of the Higgs is due to new states or whether the required level of fine-tuning in the MSSM 
is realistic, it is reasonable to consider scenarios with either a large or small splitting between the 
two stops, yielding minimal spectra in which either a single stop, or both stops and a sbottom, are 
light. 

Recently, there have been many theoretical studies on the collider signatures of light directly- 
produced third-generation squarks in the i?-parity conserving MSSM (e.g., [1, 9-16]), feeding into 
a vast experimental program devoted to hunting for such stops and sbottoms. Much less work has 
been done (e.g., [17]) on exploring the rich variety of possible signatures of third-generation squarks 
in RPV scenarios. 

One of the reasons that RPV scenarios have been under-studied is the larger parameter space - 
RPV models add several dozen new parameters to those of the MSSM. However, by restricting 
ourselves to the light stop, a general study of the RPV phenomenology is feasible. Another reason 
that i?-parity conserving scenarios have received more attention is that i?-parity avoids unacceptably 
large proton decay rates and a plethora of other disasters (for reviews, see [18, 19]). Additionally, 
i?-parity conservation makes the lightest superpartner stable, allowing it (if neutral) to serve as 
a dark matter candidate. Still, i?-parity conservation is merely a phenomenologically motivated 
assumption. An alternative to forbidding RPV couplings completely is setting many of them to 
be very small. There is no a priori reason that a technically natural hierarchy of RPV couplings, 
which satisfy phenomenological constraints, cannot be realized in nature (in fact, various theoretical 
frameworks have been shown to give rise to such situations, e.g., [20, 21]). It is therefore important 
to experimentally address all possible signatures of the RPV MSSM, especially if ongoing searches 
exclude the naturally light stop in the i?-parity preserving MSSM. 

In this paper, we will present a classification of the RPV scenarios with a light stop, using 
simplified models in which a single trilinear RPV coupling dominates, and examine to what extent 
each model is constrained by current LHC searches. We focus on the case where one of the stops 
is somewhat lighter than the other third generation squarks, so that it dominates the production, 
since typically this would be the most conservative scenario. Production of the second stop and 
a sbottom will be studied as well, although with less generality. Many of the signatures we will 
study are sufficiently simple that one could also imagine them arising outside the context of the 
RPV MSSM. In effect, the RPV MSSM can be viewed as a "signature generator" with the various 
possible final states viewed as simplified topologies [22] for generic models of new physics. We hope 
that our results will serve as guidance for designing new searches to cover the regions of parameter 
space where existing analyses are not sensitive. In addition, our results can be useful for guiding 
theoretical model building by elucidating the experimental status of the various RPV scenarios, 
very few of which have been explicitly analyzed by the experiments thus far. 

Section 2 provides a brief overview of i?-parity violation and the simplifying assumptions that 
will be made in the context of this study. Section 3 discusses the ways in which the stop may decay, 
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either directly or through other superpartners. In section 4, we use a large set of recent ATLAS and 
CMS searches to derive bounds on the various scenarios from direct stop pair production. We then, 
in section 5, study how the bounds are improved if the second stop and a sbottom are sufficiently 
light to be readily produced. Section 6 discusses the scenarios in which naturally light stops are not 
excluded by the existing searches and propose strategies which may better address these signatures. 

2 i?-parity violation in the MSSM 

When one lists all possible renormalizable couplings in the superpotential that respect the gauge 
symmetries of the MSSM, four types of "disastrous" terms appear,^ 

W 3 \k,,UL,EI + KjkLiQ.Dl + ]^>:i^,UtDpi + ^i,UH^ (2.1) 

These terms are disastrous because they all violate baryon and/or lepton number. Together, they 
can lead to proton decay and other processes which have very strict bounds. Therefore, these 
couplings are often assumed to vanish, which can be obtained by assuming the existence of the 
"i?-parity" symmetry. More generally, some of these couplings may be non-zero; one only needs 
to assume that the combinations of couplings that have stringent bounds are sufficiently small. 
Additionally, soft SUSY-violating RPV couplings are possible: 

C 3 ]^Ai,kUL,il + A^jkLiQA + + BiUK + mi;Uh\ + h.c. (2.2) 

See [19] for a thorough review of -R-parity violating operators. 

In this work, we will focus on the trilinear couplings A, A' and A", which give rise to vertices 
coupling a sfermion to two Standard Model fermions. The bilinear /ij terms in the superpotential, 
which mix the lepton and Higgs superfields, can be rotated away, adding contributions to the A and 
A' terms (primarily Aj33 and A'j33). The soft bilinear couplings, Bi and m^^, may introduce further 
mixings, but we will assume them to be negligible.^ These soft couplings can be absent at the 
messenger scale if, for example, SUSY breaking is mediated by gauge interactions. The couplings 
remain small at the weak scale if the messenger scale is sufficiently low. 

The RPV "y4-terms" of eq. (2.2), A, A' and A", will also be ignored. As these couple a stop 
to two different supersymmetric partners, the four-body decay (as the stop is assumed to be the 
lightest colored superpartner) mediated through an RPV A-term would require at least three im- 
parity violating couplings - the A-term and two terms from eq. (2.1). As the RPV is required to 
be somewhat small to be consistent with indirect experimental constraints, it is likely for these 
branching fractions to be negligible or for the decay to be displaced. If the sleptons are lighter than 
the stop, three-body decays with only two RPV couplings would be allowed, however, these would 
be competing with three-body, RP-conserving decays to the same sleptons. We will not address 
these scenarios in this work. 

^Here the Li are the left-handed lepton doublets, Ef the right-handed leptons, Qi the left-handed quark doublets, 
C/f and Df the right-handed quarks, and is the Higgs that couples to the up-type quarks; i,j,k = 1,2,3 are 
generation indices. While the couplings are potentially complex, the phase is not important because only the 
magnitude, |A|^, enters in the processes we study in this paper. 

^If this is not the case, two body-decays of Higgsinos to a charged lepton or neutrino, and a higgs, W, or Z 
become possible. The focus of this work is on the trilinear couplings, and exploration of this possibility is left for 
future work. 
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While one could consider many different kinds of hierarchies between the RPV couplings, or 
be guided by a specific hypothesis such as minimal flavor violation (MFV) [23, 21], we will base 
our study on simplified models in which a single RPV coupling, with a single choice of the gener- 
ation indices i,j, k,^ is assumed to dominate the phenomenology. These models are automatically 
viable from the point of view of indirect constraints, since bounds on individual couplings are quite 
weak [19, 24]: for TeV scale sfermions, with only a few exceptions, the bounds are O (0.1) or higher. 
Each simplified model directly corresponds to a small number of particular collider signatures, mak- 
ing the connection with experimental searches straightforward. To obtain bounds on a scenario in 
which decays through more than one RPV coupling are relevant, one may rescale our cross section 
limits by two powers of the branching fraction. 

3 Overview of signatures 

As was discussed in the Introduction, the combination of naturalness arguments with the null results 
in current LHC new physics searches suggests that one should consider spectra in which the two 
stops and at least one sbottom are much lighter than most of the other colored superpartners. This 
motivates us to make the conservative assumption that the dominant SUSY production mechanism 
is the pair production'* of the lighter stop (ti, which we will denote simply as t). In section 5, we 
will include the production of the heavier stop (^2) and the sbottom (61) and show that even if these 
particles are not very much heavier than ti, the improvements in the limits need not be drastic. 
For that reason, the limits we derive for tit^ production are realistic despite being conservative. We 
will start by considering all possible decay topologies of the stop in the RPV MSSM. 

3.1 Stop decay modes 

The stop may decay directly via an RPV coupling as 

i ij (LQD with \'-^^) (3.3) 

or 

i^jj (UDD with A'3',.,) (3.4) 

where j is a quark or antiquark (a jet, possibly a 6-jet) and £ is a charged lepton (e, /i or r). 

Alternatively, if the Higgsino or one of the gauginos is lighter than the stop and/or the RPV 
couplings that would allow direct decays happen to be sufficiently small, the stop decay may proceed 
via other (on-shell or off-shell) superpartners. The possibilities are: 

i ^ x't, x^b, ~gt, bW+, bH+ (3.5) 

with the superpartners decaying as 

X° ^ Ui^ (LIE) or ijj, ujj (LQD) or jjj (UDD) (3.6) 

X+ ^ m, ivv (LIE) or V33 (LQD) or jjj (UDD) (3.7) 

•^Barring the antisymmetry conditions A^fc — —\jik and X'l^j. = —X'^f^j, the coupUngs are independent. 
''If the UDD couphng is sufhciently large, single stop production, ds — >■ i, can also be observable. We will 
not consider it here. See [25] for a recent study of this possibility. 
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Figure 1: Left: Stop decays through an off-shell sbottom. Since the couples only to 6^ while the 
RPV couplings (apart from those allowing the stop to decay directly) involve the sbottom needs to be 
mixed. Right: Four-body cascade decays of the stop, where the intermediate particles {X = x^, or 
g, and f = D, £, or q) may or may not be on-shell. For X much heavier than the stop, diagrams with a 
helicity flip on the X propagator will dominate, as discussed in appendix A. 



~g ^ ^JJ, i^JJ (LQD) or jjj (UDD) (3.8) 

b ^ ej, uj (LQD with A^3) or jj (UDD with A'/.g) (3.9) 

Here, the sbottom decays directly through an RPV coupling (see figure 1, left), while the "inos" 
(Higgsinos and gauginos), which are not present in any of the trilinear RPV terms, decay through 
a sfermion, as illustrated in figure 1, right. As an example, with the X'123 coupling (that is, first 
generation lepton doublet, second generation quark doublet and a right handed b), one finds the 
six competing final states of t — ?• tW^* — ?• te~cb, te~^cb, tVesb and Wesb for neutral wino mediated 
decays and t — )■ bW^* buecb and be'^sb for charged wino mediated decays. 

If the stop decays most of the time into an on-shell neutralino or chargino while the RPV 
couplings involving the stop are non- negligible, these inos may preferentially decay back through 
the original stop. If only the neutralino is present (the predominantly bino case), this allows for 
four-top final states 

i^x't, t£j (LQD with A',3,) or tjj (UDD with A'3'^.,) (3.10) 

which would sometimes be easier to discover than the direct RPV decays through the same cou- 
plings. Analogously, the decay to a chargino in the wino or Higgsino case would give rise to 

i^X^b, x^ ^ ibj (LQD with A^gJ or bjj (UDD with A'3'^.,) (3.11) 

which yields events with a large number of 6-jets. Of course, diagrams where the ino decays through 
other sfermions may also contribute, with final states as in (3.10)-(3.11), but, in the LQD case, 
with the additional final states of 

X° ^ buj, x+ ^ ^tj (LQD with A^gJ (3.12) 

Multistage cascade decays that involve several inos are also possible. To keep the parameter 
space manageable, we will only treat the possibilities that may appear naturally within the wino 
system, which contains both a chargino and a neutralino, or the Higgsino system, which contains 
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a chargino and two neutralinos. If there are no light sfermions to which the inos can decay, it 
is plausible for a heavier member of the ino multiplet to transition into a lighter one (emitting a 
couple of soft SM particles - the decay products of an off-shell W, Z or h), which will only then 
undergo an RPV decay through an off-shell sfermion. If the transition is between two neutralinos, 
the RPV decay products remain the same, so there is no need to consider such cases separately. 
On the other hand, if the transition is between a chargino and a neutralino,^ 

xf-^XiW^*^^ (3.13) 

the resulting final state of the stop decay is different. Whether the chargino decays as in (3.13) 
or via a sfermion and its RPV coupling as in (3.7), is model-dependent, so we will study both 
possibilities. Similar transitions between on-shell sleptons are also possible, but we will not consider 
the possibility that they dominate over the two-body RPV decays of the sleptons. 

Even with only a single RPV coupling, several decay paths (through different mediators) can 
potentially contribute, each leading to a different final state. Which of these dominates depends on 
the masses and couplings of the intermediate particles and the available phase space. We stress that 
the i?-parity conserving decay of the stop into an on-shell ino, which subsequently decays via RPV, 
is a very realistic scenario. Naturalness suggests the Higgsinos should not be much heavier than 
~ 200 GeV (see, e.g., [1]),^ and nothing forbids the bino and winos from being light. In order to 
span the various possibilities, we will construct a separate simplified model for each RPV coupling 
and each type of mediator, with several benchmark points for the mediator's mass (both lighter 
and heavier than the stop), as will be described in more detail in the next section. 

3.2 Lifetime 

Depending on the details of the model, the decay length of the stop (or the superpartners it decays 
into) can vary from prompt (even in the case of four-body decays through off-shell mediators) to 
detector-stable. For example, a generic expression for the four-body decay rates of the stop would 
be 

3072^1;^; [^j) ^'-'^^ 

where A is the relevant RPV coupling, and Anai = gf for gauginos and the product of the relevant 
effective Yukawa couplings for Higgsinos. The resulting decay length can be written as 

. (a X 10-= ™n) (-) (-) (^) (^) (— ) (3.15) 

If the intermediate particles are not much heavier than the stop and the RPV coupling is relatively 
large, the stop will decay promptly, while for a small RPV coupling and/or heavy mediators the 
stop can easily be long-lived. 

Detector-stable stops have already been excluded by LHC searches in the natural range of stop 
masses [26, 27]. The same searches could likely set very strong limits on scenarios where the stop 
decays promptly to a detector-stable charged particle, such as a charged slepton or a chargino; 

^We assume the chargino to be heavier than the ncutrahno. Scenarios with the opposite hierarchy are possible, 
but rare. 

^In extensions of the MSSM, the Higgsinos may be heavier (e.g., up to ^ 350 GeV [8]) and stih be natural. 
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it would be useful for these searches to recast their analysis to cover such scenarios. Cases with 
a long-lived neutral particle are covered by ^^-based searches for i?-parity conserving SUSY with 
prompt decays. The intermediate case of displaced decays within the detector is less straightforward. 
Interpreting the results of the existing displaced decay searches [28-31] in the context of the variety 
of decays that we will consider here is beyond the scope of this work. Doing so reliably would require 
a close familiarity with the details of the detector and the subtleties of the particular analyses. Also, 
it is not obvious that gluinos and/or all the squarks are excluded from being as light as the stops 
in such case.^ We will therefore restrict ourselves to the case of prompt decays. We hope that our 
discussion of the various stop decay signatures will motivate the experiments to include them as 
additional benchmark models in searches for displaced decays. 

The possible longevity of the stop in RPV scenarios may allow processes through which a stop 
oscillates to an antistop (or vice versa) after hadronization [33]. As a result, the relative charges 
of the decay products of the two stops in the event can be different from the naive expectation. 
The minimum lifetime required for oscillations to occur is model-dependent, and may allow prompt 
decays. We will study scenarios with oscillations in section 4.3. 

For sufficiently light stops, there is an additional signature, the annihilation signal of "sto- 
ponium," a near-threshold stop-antistop bound state. The stoponium would decay primarily by 
annihilation if the stop decay rate is somewhat suppressed, and lead to a prompt 77 resonance at 
about twice the stop mass (see, e.g., [34-38]). However, even if stoponium is observed, studying 
the (displaced or prompt) decays of the stop itself will still be important for verifying its identity, 
the presence or absence of i?-parity violation, and potentially other properties of the model. 

4 Limits on simplified models with ii production 

As discussed in the previous section, even for a fixed RPV coupling, the stop may decay in several 
different ways, either directly to a two-body final state or through an on- or off-shell Higgsino H 
(which includes two neutralinos and a chargino), bino B (neutralino), wino W (a neutralino and a 
chargino), gluino g or sbottom b. To cover a broad range of possibilities, we construct a separate 
simplified model for each of these decay channels. For on-shell H and W, we present separate 
models for the situations in which a chargino (x^) transitions to a neutralino (x^) before decaying. 
Our simplified models for stop decays and the final states obtained in each case are listed in tables 1 
(LLE couplings), 2-3 (LQD couplings), and 4 (UDD couplings).^ These will be described in more 
detail in sections 4.1-4.2. In section 4.3, we will examine the same scenarios with stop-antistop 
oscillations. 

To determine to what extent the existing LHC searches provide coverage of these models, we 
have simulated all available 7 TeV ATLAS and CMS searches which could potentially be relevant 
and easily recast for our study.^ These are summarized in table 5. Details of the simulation are 
described in appendix B, where we also explain how we eliminate unreliable limits from searches 
probing the tails of our signal distributions. The resulting limits are presented in figures 2-5. For 
couplings related to each other by exchanging electrons and muons, or first- and second-generation 

^For a recent discussion on searches for displaced decays in the context of a specific RPV model, see [32]. 
^Throughout the text, LiQjD'^ will be referred to as LQDijfc, etc. 

^Several additional potentially relevant searches [81-86] appeared when this work was in its final stage and thus 
were not included in our analysis. 
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LLE 


mediators 


final state (of each stop) 


ijk 


first 


second 




X° ^ RPV 


121, 122 


B 
W 




Uvt 

£uvb, Uvt, (mb) 


Uvb 


131, 231, 

132, 232 


H 

B 
W 


vIIlIIr 


Urb, eriyt 

Uvt, iTVt 

tvvb, Uvt, iTvt, {Urb) 


iTvb 

Uvb, iTvb 




H 


TR 


ii/ub, irvt 


iTvb 


123 


B 


v/Il/tr 


irvt 






W 


y/h 


Tuub, irvt, {Urb) 


iTvb 




H 




iTTb, TTVt 


TTvb 


133, 233 


H 
B 


TR 

v/^l/tr 


Ivub, Tvub, Irvt, rrvt 

irvt, TTVt 


iTvb 




W 




TVUb, Irvt, TTVt, (iTTb) 


Irvb, TTvb 



Table 1: Simplified models with LLE couplings. Couplings that will be analyzed explicitly are indicated in 
bold in the first column. We denote £ = e, fi. If certain final states dominate due to phase-space suppression 
of final states with tops, they are shown in bold. Final states that are suppressed for heavy inos due to 
helicity considerations (see appendix A) are shown in parentheses. In all final states with a single charged 
lepton {£ or r), its charge is correlated with that of the stop. 

quarks (all grouped together in tables 1-4), we present results for just one representative case.^° In 
the figure captions, we indicate which of the searches from table 5 turned out to be most relevant 
for each scenario. 

4.1 Direct two-body decays 

In the absence of large hierarchies in the RPV couplings, the direct two-body decays through 
couplings involving the stop will dominate over multi-body decays that require virtual superpartners. 
In this subsection, we assume that these decays also dominate over i?-parity-conserving decays to 
lighter superpartners (if such decays are available), which may or may not be the case, depending 
on the size of the relevant couplings. The other possibility will be addressed in the next subsection. 

The A'j3^ (LQD) couplings mediate stop decays to a charged lepton (e, /i or r) and a quark 
(which may be a 6 quark). If we assume for simplicity that a single A' coupling dominates, the 
flavors of the lepton and the quark are fixed. These RPV stops are identical to leptoquarks in both 
production and decay (in the case where the branching ratio of the leptoquark to a neutrino and 
a quark vanishes). The leptoquark searches [42, 43] are directly applicable to the corresponding 
RPV stop scenarios, excluding them from the natural range of masses. However, since the third- 
generation leptoquark search [43] requires a r+b final state, signatures with a r and a light-quark 

"'^'^Exchanging electrons and muons occasionally does have an effect on the limits due to differences in the trigger, 
identification and isolation requirements and fake rates, and different fluctuations in data in analyses that have 
separate search regions for electrons and muons. However, this effect is typically small. 
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LQD 


mediators 


final state (of each stop) 


ijk 


nrst 






X° ^ RPV 


111, 112, 121, 
211, 212, 221, 


122, 
222 


g 

B 
W 


~qlylh 


^tjj, vtjj 
vbjj, itjj, vtjj, (ibjj) 


ebjj [SS], vbjj 






9 


q 


Itbj, vtbj 








H 




itbj, vtbj 


£bbj [SS], vbbj 


113, 123, 213, 


223 


B 


-qlvlh 


ttbj, vtbj 








W 


q/y/h 


vbbj, itbj, vtbj, {tbbj) 


Ibbj [SS], vbbj 






b 




ewj, vWj 








9 


q 


Ttjj, Vtjj 




311, 312, 321, 


322 


H 

B 


ql^rlTL 


rtjj, Tbjj 
Ttjj, vtjj 


Tbjj [SS] 






W 


ql^rlTL 


^bjj, Ttjj, vtjj, {Tbjj) 


Tbjj [SS], vbjj 






9 


q 


Ttbj, vtbj 








H 


bR 


Ttbj, vtbj 


Tbbj [SS], vbbj 


313, 323 




H 

B 


q/l^r/TL 


Ttbj, Tbbj 
Ttbj, vtbj 


Tbbj [SS] 






W 


ql^rlfL 


ubbj, Ttbj, vtbj, {Tbbj) 


Tbbj [SS], vbbj 






b 




TWj, vWj 





Table 2: Simplified models with LQD couplings (continued in table 3). Couplings that will be analyzed 
explicitly are indicated in bold in the first column. We denote I = e, ^; j = u, d, c, s. If certain final states 
dominate due to phase-space suppression of final states with tops, they are shown in bold. Final states that 
are suppressed for heavy inos due to helicity considerations (see appendix A) are shown in parentheses. 
For final states with a single i, r, t or W, cases with potentially same-sign dilepton events are indicated 
with [SS]. 

jet (from X'^^^ or A332) are not covered. However, searches for hadronic taus+jets+MET, and for 
the low masses the ti cross section measurements in the dilepton channel, are sensitive. These 
results are included in figure 2 (left). 

The X'^ji^ (UDD) couplings facilitate stop decays to a pair of quarks (one of which may be a 
6-quark). Searches for pair-produced particles which each decay to two jets have been done in other 
contexts by ATLAS [39, 40] and CMS [41]. In the right side of figure 2, we show our simulated 
distribution of the average reconstructed mass of the two dijets, the peak being searched for in [41]. 
We show the distribution for the stop as well as for the coloron model [88] that was actually used 
as a signal hypothesis in that CMS search. The limits derived in that search can only be translated 

^^For much of the mass range, the stop cross section is only slightly larger than our limit (see LQD332 in figure 2, 
left), which does not account for the systematic uncertainties of our simplified detector simulation. Therefore one 
should not be overly confident about our exclusion range in this case. A more dedicated experimental analysis is 
desirable. 
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LQD 


mediators 


final state (of each stop) 


ijk 


first 


second 


X > X — > V 


X° ^ RPV 
















H oiW 


i 


Ittj, ibbj 


itbj 




H ovW 




vthj 


ubbj 


131, 132, 231, 232 


W 




ittj, ibbj, utbj 


itbj, ubbj 




B 


i 


ittj 






B 


k 


vthj 






B 


q/^/h 


ittj, utbj 










ib 






H ovW 


i 


ittb, £bbb 


itbb 




H ovW 


H 


vtbb 


ubbb 


133, 233 


W 




ittb, Ibbb, utbb 


itbb, ubbb 




B 


i 


ittb 






B 


H 


vtbb 






B 




ittb, utbb 






— 





Tj 






H ovW 


i 


Tttj, rbbj 
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Tttb, rbbb 
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Ttbb 
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B 


i 


Tttb 






B 


h 


utbb 






B 




Tttb, utbb 





Table 3: Simplified models with LQD couplings through which the stop can decay directly (therefore, only 
cases with on-shell inos are considered). Same conventions as in table 2 apply. 

onto stops (after including the selection efficiencies) because these distributions have approximately 
the same shape. As we show in figure 2 (left), this search sets no limits on the stop (whose cross 

^^From our communication with CMS, we understand that future versions of their analysis may use a different 
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tbbj [SS] 




B 


q 


ttbj 





Table 4: Simplified models with UDD couplings. Couplings that will be analyzed explicitly are indicated 
in bold in the first column. We denote j = u,d,c,s. If certain final states dominate due to phase-space 
suppression of final states with tops, they are shown in bold. For final states with a single t or W, cases 
with potentially same-sign dilepton events are indicated with [SS]. The bottom part of the table includes 
couplings through which the stop can decay directly (as with the LQDz3j couplings of table 3, only cases 
with on-shell inos are considered). 

section is much smaller than that of the coloron, due to both spin and color), neither do the ATLAS 
searches at lower masses. These decays to two jets, even if one of them is a 6-jet, also receive no 
appreciable limits from any of the other searches we examined. 

4.2 Decays via other superpartners 
4.2.1 Definition of simplified models 

Let us start by considering decays, shown in figure 1 (right), that proceed through two intermediate 
particles: an "ino" X (gluino, Higgsino, bino or wino), and a sfermion / (squark, sneutrino or 
slepton). 

In addition to affecting the kinematic distributions of the decay products, the ino mass sometimes 
determines the branching ratios into the various final states. For a detailed discussion of this effect, 
see appendix A. To cover the different possibilities, we examine three cases for the mass of the 
ino: much heavier than the stop (in practice, we set rrix = 2 TeV), slightly heavier than the stop 
(mj^ = l.lmj) and on shell (with the masses discussed below). 

model for the coloron. If this changes the width of the coloron bump, then it will be impossible to do this kind of 
re-interpretation. For the ATLAS search, whose updated version [40] has appeared just recently, we did not do full 
simulation, and the presented limit assumes the cut efficiencies and width of the stop bump to be identical to those 
of the "sgluon" of [40]. Ideally, CIVIS and ATLAS would include the RPV stop as one of the benchmark models for 
which they optimize the searches in this final state. 

^•^The behavior in the transition region between light and heavy inos depends somewhat on the stop mixing angle. 
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final state 


collaboration 


£ (fb-i) 


ref. 


pairs of dijets 


ATLAS 
CMS 


0.034, 4.6 
2.2 


[39, 40] 
[41] 


leptoquark pairs 


CMS 
CMS 


5.0 
4.8 


[42] 
[43] 




ATLAS 


0.70 


[44] 


tt 


CMS 


L14 


[45] 




CMS 


2.0-2.2 


[46] 


ti + jet 


CMS 


5.0 


[47] 


ti + niT 


ATLAS 


1.04 


[48] 


{) 1 1 A /Tin"' 
I + jets + Mi!ji 


CMS 
ATLAS 


4.7 
4.7 


[49] 
[50] 


OS U + MET 


CMS 


4.98 


[51] 


OS/SS U + MET 


ATLAS 


1.04 


[52] 


SS U + MET 


ATLAS 


2.05 


[53] 


SS U 


ATLAS 


1.6, 4.7 


[54, 55] 


SS U (+ MET) 


CMS 


4.98 


[56, 57] 


SS + 6 (+ MET) 


CMS 


4.98 


[58] 


y (SS or U + 6) 


CMS 


4.9 


[59] 


3 or 4 £ 


ATLAS 


1.02 


[60, 61] 


3 ^ + MET 


ATLAS 


2.06, 4.7 


[62, 63] 


4 £ + MET 


ATLAS 


2.06 


[64] 


3 or 4 £ (+ MET) 


CMS 


4.98 


[7] 


1 or 2 T + jets + MET 


ATLAS 


2.05 


[65, 66] 


h + jets + MET 


CMS 


1.1, 4.98 


[67, 68] 


> 36 + jets + MET 


ATLAS 


4.7 


[69] 


6 (+ ^) + jets + MET 


ATLAS 


2.05 


[70] 


6 + ^ + jets + MET 


CMS 


4.96-4.98 


[71, 72] 


Z + jets + MHj i 


CMS 
ATLAS 


4.98 
2.05 


[73] 
[74] 


2-6 jets + MET 


ATLAS 


4.7 


[75] 


6-9 jets + MET 


ATLAS 


4.7 


[76] 


jets + MET 


CMS 


1.1, 4.98 


[77, 78] 


(6)-jets with 


CMS 


1.14, 4.98 


[79, 80] 



Table 5: 7 TeV LHC searches used for inferring limits. 
For the gluino, we will only treat the heavy case, since a light gluino is disfavored experimentally, 



In our models we assume the stop to be an equal mixture of left and right. 



12 




10 



,LQD332 
LQDm 
LQD232 



tlt{ 




0.10 

> 

O 0.08 
>^ 0.06 



0.04 



0.02 



100 200 300 400 500 600 700 800 
Wstop (GeV) 



0.00 



stop (UDD312) 
coloron 




200 250 300 350 400 
average mass (GeV) 



Figure 2: Left: Limits on two-body decays of the stops. The thick black curve is the stop pair-production 
cross section [87]. For decays to a lepton and a quark (LQD) we present the hmit from the CMS search 
for second generation leptoquarks [42] for the LQD232 case. A similar limit applies to the LQD131 case 
based on the search for first generation leptoquarks [42]. These limits cover the analogous cases with 
6-jets (LQD233 and LQD133) as well. We also present the limit from the CMS search for third generation 
leptoquarks [43] relevant to the LQD333 case. In the LQD332 case, which does not have a dedicated search, 
we obtain a limit from the ATLAS search for 2T-f-jets-|-MET [66] and the CMS search for opposite-sign 
(OS) dileptons-l-MET (with r's) [51], and at low masses from the it cross section measurements in the 
dilepton channel [44, 45]. For decays to pairs of jets (UDD), the limits from the ATLAS searches [39, 40] 
and the CMS search [41] are shown. None of the other searches we examined has appreciable sensitivity to 
these UDD3j/c decays. Right: in the context of the CMS paired dijets search [41], distributions (normalized 
to 1) of the average mass for the coloron model [88] used in [41] and our UDD312 stop signal, both for 
m = 320 GeV. 



and the LQD and UDD couplings, as the gluino does not couple to sleptons. As the Higgsinos 
couple most strongly to the third generation sfermions, we will only consider decays via RPV 
couplings that contain third generation indices. For on-shell Higgsinos or winos, we will assume 
rrij^ = mi — 100 GeV, which allows t — )■ hx^ but forbids t — )■ tx^ . This is a sensible choice of 
simplified model since the latter final state would be phase-space suppressed even if allowed. In the 
bino case, where the chargino is absent, we will allow t — )■ W^hx^ , with, again, m^^ = m^— 100 GeV. 
For couplings that allow the stop to decay directly (LQDzSfc and UDDSjTc), we set the bino mass 
to mj5^ = rrif: — 200 GeV, allowing for both the top and the bino to be on-shell; this allows t — )■ tx^ 
to dominate over the two-body RPV decay for reasonable parameter choices. 

In some cases, there are several contributing sfermion mediators, /, whose masses are generically 
unrelated. In particular, the masses of the right-handed sleptons, the left-handed sleptons (and 
sneutrinos), and the squarks, are independent, so their relative contributions may vary. The relative 
contributions may depend also on tan/3 (in if-mediated scenarios) since it affects the couplings of 
the different sfermions to the Higgsino differently. However, apart from kinematics, the choice of 
the dominant sfermion does not affect the possible final states or their branching ratios (so we will 
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assume all sfermions to have the same mass and set tan /3 = 40). The LLE233 and LQD323 H- 
mediated scenarios are exceptions, where different choices of sfermion mediators give rise to different 
simplified models, as indicated in tables 1 and 2. 

Additionally, scenarios where an on-shell ino can decay through the stop, as well as through 
some other sfermions, can be sensitive to the sfermion masses. Let us consider the LQD couplings 
first. Note from table 3 that for any ino, the final state for decays through hi differs from those of 
decays through t. Even though we assume the sbottom to be heavier than the stop, its processes can 
still contribute significantly if the stop is largely right-handed. Models that involve both final states 
are therefore interesting. Furthermore, in the bino and wino cases, the sleptons (and sneutrinos) 
contribute the same mixture of final states as a left-handed stop and sbottom combined (while in 
the Higgsino case, the decays through sleptons give the same final state as the decays through the 
stop). To represent the various possibilities, several cases are included in figure 5. The "mixed" 
scenarios we present were generated by forcing the winos (or binos) to decay through the left- 
handed sleptons v.Il (which is also what would be contributed by degenerate ii and 6l). At 
a typical point in parameter space, the final state composition will be somewhere between the 
stop-only and the mixed cases. It is less typical, although possible, to have the sbottom-mediated 
processes dominate. In the UDD cases (bottom half of table 4 and figure 6), for binos we do not 
separate the contribution of the stop since the final state is the same as coming from the other 
squarks. In the UDD323 Higgsino case, we present results only for decays via i since the final state 
contributed by is highly phase-space suppressed, except in the case that the chargino transitions 
to a neutralino before decaying, in which case the final state is simply the same as in the t case. 

For the purpose of event generation, we assume the sfermions / to be heavy (in practice, we 
set = 2 TeV), except when the sfermion is the stop itself, in which case we use the actual stop 
mass. The precise value of the sfermion mass does not affect the branching ratios to the various 
possible final states, and even though the details of the kinematics would have a small effect on the 
efficiencies of the searches involved, our results would still approximately apply for scenarios with 
lighter sfermions as well (which, in some cases, is actually necessary in order for the decay to be 
prompt). 

For couplings involving the right-handed bottom superfields, i.e., LQDij3 and UDDij3 (except 
for LQDi33 and UDD3j3 which allow the stop to decay directly), there are viable three-body decay 
paths for the stop through a W boson and an off-shell sbottom (see figure 1, left).^^ For regions 
of parameter space with a sufficiently mixed sbottom, this decay can dominate over the four-body 
ino- mediated decays. One may also consider an analogous diagram with a charged Higgs, 
instead of the . For simplicity, we assume that the charged Higgs is too massive for this to 
be relevant. Even if it were sufficiently light to appear on-shell, the corresponding process would 
generally be subdominant to the process (unless the sbottom mixing angle is very small), as 
the more massive is suppressed by the available phase space. 

^^For binos, we do not present the contribution of separately since the corresponding final state is the same as 
the one dominating in the mixed case. 

^^To generate these decays, we split the sbottom masses, placing them at f .5 and 2.0 TeV for the left- and right- 
dominated respectively. However, the choice of sbottom masses has no effect on the final state and a minimal effect 
on the decay kinematics. 
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Figure 3: Limits on stops decaying via other superpartners (see legend) in the presence of LLE operators. 
The first intermediate particle is either very heavy (dotted curves), 10% heavier than the stop (dashed) or 
100 GeV lighter than the stop (thin solid: x^,x^ — )• RPV, thick sohd: ^ ^ RPV). In most cases, 
the best limits are set by multilepton or same-sign (SS) dilepton searches, in particular [7, 55, 56, 59, 64]. 
Searches for OS dileptons+jets+MET [51, 52] have comparable sensitivity in the 122 and 231 heavy W 
and 123 heavy H cases, and provide the best limits on on- shell (with x^ RPV) or hght off-shell 123 H 
cases. In the 233 H-fji cases with x^jX^ ~^ RPV, the limits are set by the ti cross section measurements 
in the dilepton channel [44-46] (for heavy together with [56, 59]), with OS dileptons [51, 52] and single 
^+jets+MET [49, 70] also being important. Searches for 6-jets+MET [68] and 1 or 2r+jets+MET [65, 66, 
51] are comparable to the SS dilepton search in the 123 heavy W case and provide the best limits on the 
233 heavy W case. 



4.2.2 Discussion of the results 

The hmits we obtain for the various scenarios are presented in figures 3-6. The different mediator 
assumptions are represented by different colors, as indicated in the legends. The style of the curve 
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Figure 4: Limits on stops decaying via other superpartners (see legend) in the presence of LQD operators 
(continued in figure 5). The first intermediate particle is either very heavy (dotted curves), 10% heavier 
than the stop (dashed) or 100 GeV lighter than the stop (thin solid: x^,x^ — )• RPV, thick solid: — ^ 
X^ —7- RPV). For decays via g, B and H (except for cases mentioned below), the limits are set by 
searches for SS dileptons [56, 58] and h' [59], except in 321 and 323 cases at very low rristop where the most 
sensitive analyses are 6-jets+MET [68] and tt cross section measurements [44, 45]. In H-Ur/^L cases with 
X*^) ~^ RPV, the most sensitive analysis is the ti cross section measurement in the I + Th channel [46]. 
In the 323 H-bji case with x^ — )• transitions, the limit is set by searches for 6+jets+MET [68, 69]. 
In most of the W cases, the best limits are set by the (5-)jets+MET searches [75, 68, 69]. In 221 and 

RPV, ^(+6)+jets+MET [49, 70], and «-like searches [47, 48] are comparably 



123 cases with x ,X 

important (except for heavy W). In 221 and 123 cases with x^ — )■ x^ transitions, the limits are set by SS 
dileptons [56, 58, 59]. For 6- mediated scenarios, the best limits are set by SS dileptons [56] and (in the 123 
case) multileptons [7, 63], and (in the 323 case, at low rn-stop) ti cross section [44, 45]. 
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Figure 5: Limits on stops decaying via on-shell inos (see legend) in the presence of LQD operators through 
which the stop can decay directly. The Higgsinos and winos are assumed to be 100 GeV lighter than the 

RPV) and the binos 200 GeV hghter than the 



stop (thin curves: jciX ~^ RPV, thick curves: x 
stop. See text for details about the meaning and interpretation of the t, hi and mixed cases. For H and W 
mediators with x^, or B mediators, the best limits on the t-dominated cases are set by the searches 

for SS dileptons [56, 58] and b' [59], while for the other sfermion choices presented, most powerful are the 
searches for 6-jets+MET [68, 69]. For H or W mediators with — )• RPV, the limits on the ^-dominated 
cases (and for 232 and 233 also the mixed cases) are set by the ti cross section measurement in the dilepton 
channel (without taus [45] for 232 and 233, and with a tau [46] for 332 and 333); in the 333 case the search 
for 6-jets+^+MET [72] is close to extending the limits to higher masses. In the 6i:,-dominated cases of 
H/W with RPV, the searches for 6-jets(+^)+MET [69, 72] are closest to setting Umits. 



describes the mass of the mediator: dotted for a mediator much heavier than the stop (which is the 
only option considered for a gluino and a sbottom), dashed for a mediator 10% heavier than the 
stop, and sohd for an on-shell mediator. For Higgsinos and winos, the on-shell mediator scenario is 
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Figure 6: Limits on stops decaying via other superpartners (see legend) in the presence of UDD operators. 
For 212 and 213 couphngs, the first intermediate particle is either very heavy (dotted curves), 10% heavier 
than the stop (dashed) or 100 GeV lighter than the stop (thin solid: x^,x^ — )• RPV, thick solid: — >■ 
X^ — )• RPV). For 312 and 323 couplings (through which the stop can decay directly), the Higgsino (bino) 
is taken to be 100 GeV (200 GeV) lighter than the stop. In the 212 and 213 cases, the search for b- 
jets+MET [68] was close to setting limits on decays via an off-shell g oi B, or any H with x^,x^ — )• RPV, 
while for 5-mediated decays, the limit is set by the ti cross section measurements [44, 45]. For the 312 and 
323 couplings, the limits on the -ff- mediated cases with x^ ^ x^, and S-mediated cases, are set by SS 
dileptons [56, 58, 59]. 



further separated to cases where ^ {RPV} (thin solid line) and where, via the process (3.13), 
X+ — )■ x'^ — )■ {RPV} (thick solid line). In the latter case, we conservatively assumed the chargino 
to be only 5 GeV heavier than the neutralino so that the decay products of the off-shell W are 
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essentially undetectable.^^ For decays through on-shell H or W, we do not present limits for 
mstop < 200 GeV since charginos lighter than 100 GeV are excluded by LEP;^^ in on-shell B cases 
we include neutralinos as light as 50 GeV. For 4-body decays through off-shell B and g mediators, 
we do not present limits for mstop < 200 GeV since the final state always contains a top, and with 
an off-shell top the five-body decay would likely be displaced (assuming in the B case that the 
sleptons are not much lighter than 200 GeV). 

For LLE couplings (table 1, figure 3), in almost all cases, the lower bounds on the stop mass are 
as high as 600-700 GeV, beyond the natural range for stops. These strong limits indicate that the 
experimental coverage of these scenarios is very good, as even before taking into account branching 
ratios, acceptances, identification efficiencies and cuts, a 750 GeV stop would have only ~ 10 events 
in 5 fb~^ of 7 TeV data. The limits are so strong because all of the events contain either four or 
more charged leptons, or at least two charged leptons and several neutrinos (a source of I^t)- The 
weakest limits for each coupling are obtained when the mediator is a heavy wino, or a Higgsino 
of any mass decaying via tr (except for cases with — )■ x° transitions), where the dominant 
signature is a pair of opposite-sign leptons, two 6-jets and I^t (in appendix A, we explain why this 
final state dominates for winos much heavier than the stop). 

For LQD couplings that do not involve the stop (table 2, figure 4), the limits are weaker than in 
the LLE case because there are fewer leptons. For light-lepton LQD couplings (221 and 123), gluino, 
bino, Higgsino, and sbottom-mediated decays, as well as scenarios with on-shell — )■ x° transitions, 
have relatively strong limits from SS dilepton searches. Indeed, decays of neutrahnos and gluinos 
can produce leptons of either sign. Decays via the sbottom or the chargino will produce opposite- 
sign leptons, but additional leptons can arise from leptonically-decaying Ws (in the sbottom case) 
or tops (in a chargino case where a top is produced in the chargino decay, as happens in the 
123 Higgsino scenario). On the other hand, SS dilepton events are rare in the wino scenarios since 
decays via x^ + b dominate over decays via x^ + t (due to phase space). Furthermore, hard neutrinos 
(i.e., ]^t) are available only at the price of losing a lepton. Consequently, the limits on the wino 
scenarios are relatively weak (with the exception of the x^ — )■ x° scenarios mentioned above). 
The picture changes drastically for couplings where the lepton is a tau (321, 323). Hadronically 
decaying taus have low identification efficiencies and high fake rates, while leptonically decaying 
taus are suppressed by their branching ratios. As a result, the signal gets spread over many different 
final states, making the limits on all the scenarios weak. 

For LQD couplings that involve the stop (table 3, figure 5), the inos have multiple decay paths 
with branching fractions that depend on the mediating sfermion masses, as discussed in the previous 
subsection. The decays of the bino through the stop give rise to events with two leptons and four 
tops, leading to very strong limits from SS dileptons. However, if the bino also may decay via other 
superpartners, the decays ^ ^bj dominate over x° Hj due to phase space suppression from 
the heavy top. This leads to final states with only a single ti pair and no leptons, resulting in very 
weak limits. In Higgsino and wino scenarios with X° {R-PV}, events with two leptons 

and two tops, from the decays of x° through the stop, again allow strong limits to be set by SS 
dileptons. The limits remain strong when the neutralino can decay through other superpartners 
due to the searches for multiple 6-jets-|-MET, unlike in the bino case where these searches are not 
sufficiently powerful. In Higgsino and wino scenarios with x^ — {RPV} decays, the dominant final 

^®For splittings > 20 GeV, transitions (between the chargino and neutrahno, or between the two neutrahnos of the 
H case) wiU sometimes add detectable objects to the events, and in some cases this will make the discovery easier. 

^^However, models with stop-chargino splitting smaller than our benchmark value of 100 GeV are possible, so 
scenarios with even lighter stops are also possible. 
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Figure 7: Limit on t — )■ T~^j decays (LQD332) in the presence of stop-antistop oscillations. The limit is 
set, with approximately equal strengths, by the searches for SS dileptons [56] and 2r+jets+MET [66]. 



states contain either two opposite-sign leptons (and jets, but no MET) or a ti pair (accompanied 
by jets and MET), so the limits are weak. 

For UDD couplings (table 4, figure 6), the limits on all the 212 and 213 cases (and the 312 and 
323 Higgsino cases with x"^ — )■ x°) are very weak. Indeed, the signatures tt+jets, Vr"'"iy~-|-jets, 
or just jets, have large SM backgrounds, and there are no searches that make use of the large jet 
multiplicities or the extra 6-tags that are available in these 1ow-^t stop events. On the other hand, 
for 312 or 323 couplings, the 4-top events of bino-mediated decays, and the same-sign di-top events 
of Higgsino- mediated decays with x"^ — )■ {RPV}, allow for limits to be set by SS dilepton searches. 

The least constrained cases, and possible methods to search for some of them, will be discussed 
in more detail in section 6. 

4.3 Scenarios with stop-antistop oscillations 

Sufficiently long-lived stops (i.e., F^ < Aqcd) will hadronize. If the stop hadronizes into a neutral 
"mesino," it may oscillate to an antistop before decaying [33]. The oscillation period is determined 
by the rate of the relevant fiavor violating processes, which is model-dependent, however it is 
plausible to have a situation in which the stop oscillates and still decays promptly. Oscillation 
can enhance the detectability of the stops, due to a higher abundance of same-sign lepton pairs in 
cases where the leptons would have opposite signs otherwise. Our implicit assumption so far has 
been of a regime where the oscillations are not important. Now we will examine the situation with 
oscillations. 

The probability for a stop to form a neutral mesino is | < /o ^ ^ (the uncertainty is due to the 
strange quark) [33] . The probability for a neutral mesino to be found in the anti- mesino eigenstate 
at the time of the decay reaches P = | in the rapid oscillation regime. To illustrate the possible 
effects of oscillations, we will assume the maximal possible value for Pqsc = 'Pfo = \- This implies 
there is a 3/8 probability for an event to oscillate to tt or t*t*, allowing for more readily produced 
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Figure 8: Limits on LQD scenarios with stop-antistop oscillations. In most cases, the best limits are set 
by SS dilepton searches [55-59]. However, the W scenarios (except for light W in the 221 and 123 cases) 
are still constrained primarily by (6-)jets+MET searches [68, 69, 75] as in the case without oscillations. 



SS dileptons. 

The resulting limit on the least constrained two-body LQD decay, t — t- r^j, is shown in figure 7. 
Except at low masses, it has not improved relative to the case without oscillations (figure 2) - the 
limit from SS dileptons is comparable in strength to the limit from the search for 2r+jets+MET. 
The limits on the two-body UDD decays (not shown) do not improve either as their final states do 
not contain leptons. 

The limits on decays through the inos or sbottom are shown in figures 8 and 9, where we have not 
included LLE scenarios since the limits there are usually very strong even without oscillations, nor 
have we included the cases in which the stop decays to an on-shell ino, assuming that the stop would 
be too short-lived to hadronize and oscillate. As comparison with figures 4 and 6 demonstrates, SS 
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dilepton searches improve limits significantly in many of the cases. In fact, these searches now set 
the best limits in all cases, except for LQD heavy-wino mediated decays (where final states with 
leptons are suppressed by helicity considerations) and LQD321/323 light-wino mediated decays 
(where 6-jets+MET final states yield stronger limits than final states with SS taus). 

5 Limits on simplified models with ti, bi production 

Having derived limits under the conservative assumption that they are dominated by tit^ pro- 
duction, we now examine scenarios containing the additional production of )f:2^2 bib^^^ (since 
naturalness also predicts t2 and bi to be light). Let us consider a case where the masses of t2 and 
bi are 

rrii^ = = mi^ + 100 GeV (5.16) 

and they decay as 

ta^tiZ, bi-^iiW- (5.17) 

This leads to the same final states as found in the tit\ case, but with additional ZZ or W^W~ . 

As shown in figure 10, UDDSjA; scenarios where stops decay to pairs of jets, which are completely 
unconstrained with tit\ production alone (figure 2), now obtain limits from searches for multileptons 
and Z-l-jets-l-MET.^^ For LQD decays to rj, the lower bound on the stop mass is still set by the 
same searches as before, even though the searches for SS dileptons and multileptons strengthen the 
exclusion at lower masses. 

^®If the sbottom has a significant right-handed component and the UDD coupling A'/23 is relatively large, single 
sbottom production may be observable, as has been studied in [89]. We will not explore this possibility here. 
^^The utility of the CMS multileptons search for scenarios of this type has been also studied earlier in [17]. 
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Figure 10: Limits on scenarios with ^2 and bi (both 100 GeV heavier than ii) with two-body decays 
of ti — compare to the ti-only limits of figure 2. In the UDD cases, the limits are set by the searches 
for multileptons [7] and Z+jets+MET [73]. In the LQD case, the limits are set by the searches for 
2r+jets+MET [66] and OS dileptons+MET [51], as with ti-only production, except at low masses where 
they are strengthened by searches for SS dileptons [56] and multileptons [7] . 



The results for decays via off-shell intermediate particles are shown in figures 11 and 12.^*^ The 
limits on scenarios that were least constrained with ti alone (LQD 321 and 323, and all the UDD 
cases) improve the most (compare to figures 4 and 6), in some cases even to the extent of almost 
excluding the stops from the natural mass range. 

However, in the presence of inos lighter than t2 and bi, the decays (5.17) can be highly suppressed 
due to the possibility to decay directly to the inos via 

k ^ bx+, (5.18) 

^1 -> bx\ tr (5.19) 

Such decays of t2, especially to inos that are lighter than ti, would typically lead to signatures similar 
to those arising from the decays of ti via the same ino, so they would effectively just increase the size 
of the signals we studied in section 4 by a small amount. For the sbottom, the dominant decay would 
typically be hi — )■ leading to a signature similar to that of ti — )■ bx^ followed by a transition of 

X~^ to as in (3.13). We have already included simplified models with the corresponding final states 
in section 4.^^ While studying these classes of scenarios in full generality would involve scanning 
a multidimensional parameter space, it is clear that the signatures and the corresponding limits 
would often be similar to those of the ti-only scenarios from section 4. Therefore, there is a stronger 
motivation for designing searches for the more conservative ti-only scenarios. This is relevant also 

^°LLE limits are not presented as most of them are very strong even in the case with only a single stop. 

^-'^ However, in the Higgsino case, the decay bi -> tx^ may dominate (to the extent allowed by phase space) due to 
the large Yukawa coupling of the top. 

^^In the bino case, which does not have a x^, the analogy does not work. However, bino decays almost always 
have the same final states as wino and/or Higgsino decays with the same RPV coupling. 
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Figure 11: Limits on scenarios with t2 and hi (both 100 GeV heavier than ti) in the presence of LQD 
operators — compare to the fi-only case of figure 4. In all the g-, H-, B-, and 6-mediated scenarios, as 
well as 221 and 123 scenarios with light W, the strongest limits are obtained primarily from searches 
for SS dileptons [56, 58] and b' [59]. In the other VF- mediated scenarios, the limits are set mostly by 
(5)-jets+MET searches [68, 69, 75], similarly to the ti-only case, with comparable limits obtained from 
searches for Z+jets+MET [73, 74] at low mgtop- 



to scenarios with off-shell mediators (the limits on which were presented in figures 11-12) since it 
is quite plausible for the t2 and bi to be much more massive than in (5.16). With a significantly 
reduced production cross-section, their final states can be irrelevant relative to those of ti. 

One should also consider scenarios in which two-body RPV decays of the sbottom dominate 
over its decays to hghter superpartners. Sbottom decays via LQDiSS operators result in a ub 
final state, which is covered by searches for bi — )■ 6xi with a stable and massless Xi [90] o^' third 
generation leptoquarks [91], where the current lower bound on the sbottom mass is 490 GeV [90]. 
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Figure 12: Limits on scenarios with t2 and bi (both 100 GeV heavier than ti) in the presence of UDD 
operators — compare to the ti-only cases of figure 6. The limits are set by the search for SS dileptons [56] 
and (except for the 6- mediated case) the b' search [59] and (in the 213 case) the search for SS dileptons 
with 6-tags [58]. 



If the sbottom has a sizable right-handed component, decays to it or rt will also be present and 
are likely to be detectable in SS dilepton or multi-lepton searches. In the LQDzS/c {k ^ 3) cases, 
the sbottom decays to uj. Even though it is a leptoquark-like final state, it is not covered by the 
existing leptoquark searches since those require at least one of the leptoquarks to give a charged 
lepton [42, 92, 93], or the jets to be 6-jets [91]. If the sbottom is partially right-handed, its decays 
through UDDij3 {i ^ 3) would give dijet pairs, a final state we have already analyzed in section 4.1. 
On the other hand, UDD3j3 would give rise to tops (pairs of tj resonances), joining the tt-f-jets 
signatures of some of the stop decays from table 4. In the LQDzj3 (i,j 7^ 3) case, the sbottom 
has and z/j final states, a signature covered by the standard leptoquark searches [42, 92, 93]. 
However, the combination of rj and z/j final states in the LQD3j3 (j 7^ 3) case in not covered since 
the third generation leptoquark searches assume the jet to be a 6-jet [43, 91]. Finally, one should 
note that if the sbottom is lighter than the stop, the stop may decay as ti — )■ h^^'^*\ If the h\ 
then undergoes a direct RPV decay, the overall stop decay final state will be similar to the decays 
through an off-shell sbottom from section 4. 

6 Summary and discussion 

Motivated by naturalness of the MSSM Higgs potential, which favors light stops, we have spanned 
the i?-parity violating stop decay topologies. We considered most plausible prompt^'^ decays through 
the various possible mediators in the MSSM. We utilized the ansatz of single RPV coupling dom- 
inance to pose the RPV parameter space in a way amenable to a comprehensive, signature-based 

^■^It is important to remember that the decays of the stop may be displaced, since the RPV couplings may be very 
small. In such a case, our limits do not apply, though our classification of the various possible decay signatures may 
be useful for ensuring that future searches for displaced decays take these possibilities into account. 
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exploration. Considering both the direct pair production of a single stop and of two stops and a 
sbottom, we recast existing experimental searches in the 7 TeV LHC data to derive bounds on these 
scenarios. 

We found that natural stops are strongly constrained from decaying via LLE operators, moder- 
ately constrained from decaying via e or /i LQD operators, and quite unconstrained from decaying 
via LQD operators with taus or UDD operators. Including the production of the second stop 
and a sbottom or stop-antistop oscillations improves limits in some scenarios, but neither of these 
factors has to be present. Let us therefore focus on the more conservative titl scenarios without 
oscillations, and explore the cases in which a naturally light stop within the RPV MSSM is still 
allowed. In what follows, we will discuss search strategies that could address these unconstrained 
or underconstrained final states. 

The various leptoquark searches have placed very strong limits on the two-body stop decays 
through LQDiSk operators. An exception is the rj final state (where the jet is not a &-jet), which 
does not have a dedicated search despite the simplicity of the signature. Such a search would likely 
be able to set strong limits on stops decaying in this manner, since even the generic searches for 
rs in association with jets and MET are able to set moderately strong limits on this signature (see 
LQD332 in figure 2). 

Some of the most difficult final states to address experimentally are stop decays into pairs of 
jets. As we have seen (figure 2), the existing searches for paired dijet resonances [39-41] currently 
place no limits on the stop mass. One simple modification which may allow for sensitivity to the 
stop is to require one or two 6-tagged jets (i.e., to form jb resonances). Indeed, two out of the three 
UDD couplings that could facilitate such a decay include a 6-quark. As additional motivation, these 
couplings dominate if the structure of the RPV couplings is determined by MFV [21]. The presence 
of 6-jets may also be useful for triggering, allowing for sensitivity at low masses, even with high 
instantaneous luminosities. Even without 6-jets, it may be possible to save bandwidth by recording 
partial event information as has been done recently in the CMS search for low-mass (single) dijet 
resonances [94]. It may also be useful to utilize multiple triggers to capture more signal events 
into the sample, to use pre-scaled triggers, or to extend the analyses down to masses where the 
trigger efficiency is not fiat. More generally, it would likely be beneficial to optimize the searches 
for the stop, rather than the sgluon or the coloron, which differ in the size of the cross section, the 
invariant mass distribution of the pairs (i.e., the boost of the stops, on which these searches cut), 
and radiation (which differs between quarks and gluons and affects combinatoric ambiguities and 
the signal shape). 

While the signature of dijet pairs is one of the most well-known examples of difficult stop decays 
in the RPV MSSM, there are many other difficult signatures which are very well motivated. For 
example, it is quite plausible for an electroweak gaugino or the Higgsino to be lighter than the stop - 
a natural MSSM even requires the Higgsino to be fairly light. The -R-parity conserving decays of 
the stops to these particles can then easily dominate over the direct RPV decays of the stop. The 
gaugino or Higgsino would then undergo an RPV decay through a diagram with an off-shell stop 
or another sfermion. As we have seen, this and other types of scenarios lead to a diverse spectrum 
of possible final states, the current limits on many of which are very weak or non-existent. 

In table 6, the scenarios for which the limit on the stop mass does not exceed 500 GeV are 
classified according to the dominant final states for the whole event. Why do these signatures 
remain elusive? Unsurprisingly, they typically contain: multiple jets, taus, little or no I^t, and no 
more than two leptons with any sizable branching ratio. However, while some of the cases may be 
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Final state 


6-jets 


Scenario(s) 


{T+j){T-j) 





LQD332 


{n){n) 


0, 2 


UDD312/323 


8j 


4, 6 


UDD312/323 with H; UDD213 with 


£~^£~ + 6j 


2, 4, 6 


LQD232/233 with H /W decaying via i 
LQD221/123 with W 


r+T^ + 6j 


2, 4,6 


LQD332/333 with H /W decaying via t 
LQD321/323 with H-Vr/fL or W (with or without x°) 




2,4 


LQD321/323 with H-Vr/^L or W, with x° 


+ 6j 


2, 4 


UDD212/213 with g/B; UDD213 with H 


it + 4j + MET 


2, 4, 6 


LQD321/323 with g/B; LQD323 with H-Ir 
LQD232/233/332/333 with H/W/B decaying via 6^ 


{tt or tt) + 6j 


4, 6 


UDD312/323 with H° 


tt + 2t + 4j 
tt + T + Aj + MET 


2, 4 


LQD321/323 with g/B; LQD323 with H-Ir 


T+T-W+W" + 2i 
r + W+W- + 2j + MET 
W+W- + 2j + MET 





LQD323 with b 


4 tops + 4j 


4, 6 


UDD312/323 with B 


6j + MET 


2,4 


LQD221/123/321/323 with W 
LQD321/323 with W'^ 
LQD323 with ^ H° ^ Ir 


^ + 6i + MET 


2,4 


LQD221/123 with W 


T + DJ + Mrlji 




LQD321/323 with W (with or without W^) 
LQD323 with ^ H° ^ Ir 


e-^£- (or ^+r-) + 26 + MET 


2 


LLE233 with H-tr 


T+T- +2b + MET 


2 


LLE233 with H-tr; LLE123/233 with heavy W 


+ 4j 





UDD213 with b 



Table 6: Dominant final states in scenarios for which the coverage is insufficient (for rustop ^ 500 GeV). 
See tables 1-4 for more detailed descriptions of the scenarios mentioned. As before, couplings related by 
interchanging electrons and muons, or first and second generation quarks, are listed just once. The second 
column indicates the possible number of 6-jets in each scenario (including those coming from top decays, 
where relevant). It is assumed that an on-shell chargino decays via RPV (rather than transitions to a 
neutralino), except where explicitly noted otherwise. 

genuinely difficult, many of these final states have unique characteristics that are not exploited in 
the existing searches: 

• Most of the final states contain at least two 6-jets. This happens because the decays of the stop 
through a chargino (neutralino) always produce a bottom (top). In some of these cases there 
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are even four or more 6-jets overall. These scenarios provide important motivation for includ- 
ing search regions with a large number of 6-tagged jets. So far, only a small number of final 
states have been studied with three 6-tags (jets+MET [68, 69, 80], lepton+jets+MET [72], and 
SS dileptons+MET [58]) and there have been no searches requiring four or more 6-tags. Some 
final states have not been searched for with 6-tagging at all. Aside from the (quite specific) 
search for third generation leptoquarks [43] and one of the ti cross section measurements [46] , 
there have been no new physics searches requiring both hadronic taus and 6-jets, which would 
be relevant for a large fraction of the scenarios with weak or no bounds in table 6. It would 
also be interesting to explore to what extent the requirement of an unusually large number of 
6-tags can replace the I^t requirement for both triggering and background rejection purposes. 
We should note that reduction in the signal efficiency due to 6-tagging (or other requirements) 
would typically not be an issue as light stops have huge production cross sections. 

• Many of the final states contain a ti pair (where the tops may or may not be on-shell) or a 
W^W~ pair (with extra jets). Therefore, searches based on ti cross section measurements 
(without harsh cuts on jet pts, etc.) could potentially constrain some of these scenarios. 
The stop signal will not stand out on top of the uncertainty on the ti cross section unless the 
stops are lighter than ~ 150 GeV (which is still possible in some cases). However, requiring 
multiple additional jets or even extra 6-tags (which are available in many of the scenarios), or 
looking at the invariant mass distribution, is likely to be helpful. 

• Some of the stops decay into final states that mimic top decays. For example, t — > Ivvh 
or Tuub, in LLE233 scenarios with the Higgsinos decaying via f^, is similar to the leptonic 
decay of the top (t — )■ iub or rub). This, again, motivates searches based on the ti cross 
section measurements (in the dilepton channel). These can be extremely efficient since the 
ti background is suppressed by the dileptonic branching ratio while the stop signal is not. 
Indeed, we have found several cases in this study where the ti cross section measurements 
turned out to be more sensitive than any existing new physics searches that we examined 
(despite the fact that we used cross section measurements with 1-2 fb~^ while most of the 
new physics searches were based on 5 fb~^ of data). In some cases, such as the decays of on- 
shell winos or Higgsinos via the stop in LQD scenarios, the stop final states (e.g., t — )■ ibbj) 
do not even contain neutrinos, but the contribution of I^t from the multiple jets in the event 
(jet mismeasurement, neutrinos from b decays) is sufficient for passing the very mild cuts 
of the ti cross section measurements. It is clear that with further optimization of the selection 
criteria, the limits on such dileptonic tt-like scenarios, many of which are not yet sufficiently 
constrained, can be improved significantly. 

• Leptons produced in the RPV decays of stops or inos are generally hard. Designing searches 
with hard cuts on lepton ptS may be useful for improving limits. An example of this gap in 
coverage can be seen in the LQD221 scenario with on-shell winos, where the dominant process 
is t — bW~^, — )■ vjj or l^jj. The strongest bounds there come from single £+jets+^T 
searches, while a quarter of the events contain two very high px leptons (and many hard 
jets). Utilizing hard leptons in such a case may be more efficient than I^t- This could be 
effectively implemented by a cut on the variable At = ^^Pr, or simply requiring the leptons 

^^Such tt-like searches can also be useful for light stops (toj < rrit) decaying to W , b, and an invisible particle in 
i?-parity conserving scenarios [9, 5]. 
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to pass harder cuts. In fact, the combination of high leptons with very httle I^t could 
be an effective discriminator against the dileptonic ti background. Additionally, multilepton 
searches in high At regions could provide even tighter limits on some of the simplified models 
already constrained, increasing the reach for scenarios in which decays through lepton-rich 
RPV couplings (LLE) compete with less spectacular, but more common decays. 

• Final states without hard neutrinos lack significant I^t, making conventional SUSY searches 
inefficient. However, many of these final states contain one or two leptons and a large number 
of jets. Such scenarios may be accessible by replacing the large I^t requirement with the 
requirement of a large jet multiplicity, as has been studied in [95]. 

• The structure of RPV couplings can easily introduce a preference for one lepton flavor over 
another. This means that keeping separate search regions for the different lepton flavors can 
improve sensitivity. Furthermore, the lepton flavor universality of most of the SM backgrounds 
can be used for data-driven background estimates (e.g., doing a measurement in the electron 
channel for estimating the expected number of events in the tau channel). 

It is our hope that the experimental community will take these suggestions into consideration 
in order to maximize the potential for constraining or discovering the light stop of the i?-parity 
violating MSSM, or any new physics beyond the standard model for that matter, at the LHC. 
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A Dependence on the ofF-shell ino mass 

In this appendix, we discuss the effect responsible for the dependence of the branching ratios on the 
ino mass for decays via an off-shell ino. In particular, we will explain why some final states (those 
presented in parentheses in tables 1 and 2) in VT-mediated decays disappear when the wino is much 
heavier than the stop. The effect is present for other ino mediators as well, but is less dramatic. 
The essence of the effect, in terms of the right diagram in figure 1, is that the ino propagator can 
either preserve or flip the helicity, and since the helicity-flipping propagator is proportional to the 
ino mass, it dominates the process for heavy inos. 

As an example, consider the VT-mediated decays with LLE couplings (figure 3), where the effect 
on the limits is most significant. Processes with tops are phase-space suppressed for any wino mass, 
so the dominant mediator is the charged (rather than the neutral) wino, i.e., X = x"*" ~ in 
the notation of figure 1 (right). The wino can decay through either a left-handed slepton or a 
sneutrino (/ = ii or u, with ^ fny), so the final state can be either hvll^*^ or W+v^*^ (with 
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the or i> decaying further via RPV). As winos only couple to SU{2)i doublets, the h is left- 
handed. Without a helicity flip, the helicity of /i would need to be right-handed, which allows it 
to be but not v. With a helicity flip it can be i/, but not . As a result, since the helicity- 
flipping diagram is proportional to m^, final states with v dominate for irty^ ^ ^i- This results in 
r(t — 7- buui^) ^ r(t — 7- bi~^i~^i~) for heavy winos, weakening the limits relative to cases with light 
winos where the rates of the two processes are comparable. Similarly, for the LQD operators and 
heavy winos, T{t ^ bujj) S> r(t — )■ b£~^jj). 

For the sake of completeness, let us also discuss the effect for processes mediated by the neutral 
wino even though they are suppressed relative to the """-mediated ones. The helicity-fiipping 
diagram produces /i = i~ or u, while in the helicity-conserving diagram /i = £^ or V (assuming 
again f = or u)."^^ As a result, the leptons from the stop and the antistop (if both decay through 
this process) will always have opposite signs for m^y ^ m^, but will sometimes have same signs 
for lighter winos (note though, depending on the case, additional leptons may be coming from the 
sfermion decays). 

For the Higgsinos, bino, and gluino, which couple to both right- and left-handed fermions, 
one can construct helicity-fiipping diagrams with either outgoing helicity (as long as the stop is 
somewhat mixed), so the effect is typically less drastic than it is for the winos. 

B Details of simulation and limit computation 

We use the RPVMSSM model [97] of FeynRules [98] to define the RPV couplings for Mad- 
Graph 5 [99] via the UFO interface [100]. To make the generation of our (up to) 2 — )■ 10 pro- 
cesses feasible, we generate events for stop pair production (2 — i- 2) and the various stop decays 
(1 — )■ 2,3,4,5) separately and combine the resulting LHE files while taking the boosts and the 
color connections of the stops into account. The combined events are showered and hadronized 
in Pythia 8 [101, 102] and further processed with a private detector simulator (using the anti-fcr 
jet algorithm from Fast Jet [103]), which uses truth MC information and includes geometric ac- 
ceptances of the various particles, jet energy resolution (based on [104]), identification of 6-jets and 
hadronic r candidates, and computation of isolation variables for leptons. We then apply trigger 
efficiencies, lepton identification efficiencies, lepton isolation requirements and 6-tagging efficiencies 
relevant to each search, to the extent that details about them are provided in the experimental 
publications or obtained via other means. For the scenarios of section 4.3, we also apply a prob- 
ability for leptons from each stop decay to reverse their charges due to a stop-antistop oscillation. 
The resulting events are passed through the analysis cuts. We then compare the NLO-I-NLL pro- 
duction cross section [87] with the 95% CL excluded cross section (the limit on a x e divided by the 
simulated efficiency) for each search region. In cases where the limits on a x e are not provided in 
the experimental papers, we compute them with the frequentist method [105] using the provided 

■^^The consequences are perhaps surprising. At first glance, one might assume that, for example, the decays 
t — > tx'^* — > t{e^Vfj^T^ vs. e^f^T+) (for A123) have equal branching fractions because the mediator is neutral. This 
would be correct for an on-shell x° but is not generally the case for off-shell inos. For m^o 3> mi, the helicity flip 
will make te^Vf^T~ dominate over the final state in which the decay products of the neutralino are replaced with 
their antiparticles. A similar effect has been discussed in [96]. 

^^Scenarios with t2 and 61 production are generated in a similar way, with the or 6is decayed down to tis in the 
first stage, and then combined with ti decays. For scenarios with chargino-neutralino transitions, the decay down to 
the neutralino is included in the first stage, and then combined with or decays. 
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backgrounds (and their uncertainties). The single search region giving the best hmit is used in each 
case. 

We have vahdated our detector simulation and analysis code on signal models that were used 
by the experimental analyses in cases where the experiments provided their simulated events yields 
for examples of such models. The previous version of our detector simulation code has been also 
validated in [5]. Typically, our event yields agree with those quoted in the experimental papers to 
within ~ 30% (although, in a few cases the discrepancy is about a factor of 2).^^ The reader may 
shift our exclusion curves by these amounts to estimate by how much such uncertainties may be 
affecting the mass limit in each case. When specifying which searches set the best limits, we have 
listed not only the searches that turned out to be the most powerful according to our simulation, 
but also those that had comparable power within this uncertainty. 

We have not explicitly taken the systematic uncertainty on the signal efficiency into account. 
In most cases, this can indeed be neglected since the overall uncertainty is typically dominated by 
the background. However, the uncertainty on the signal efficiency becomes very important when 
the analysis cuts are such that the signal is coming from the tails of the distributions (of I^Ti H^, 
etc.). The tails are problematic because higher-order QCD corrections and/or imperfect modeling 
of the detector may change them dramatically. The multiplicity of scenarios and searches that we 
cover here does not allow us to analyze the range of validity in each case in detail (we also do not 
check whether there are any cases in which the contribution of the signal to the control regions of 
the searches has a significant effect on the background estimate). However, we have used the size 
of the signal efficiency e as a rough proxy of sensitivity to the tails (ignoring the fact that some of 
the efficiency reduction comes from branching ratios) and included only limits that are based on 
efficiencies above emin = 10~^. More specifically, we have implemented this threshold by modifying 
the excluded cross section as cr — )■ crexp(emin/e)- For scenarios that combine the ii, ^2, and bi pair 
production processes, in order to take into account the qualitatively different properties of these 
three samples and their unequal cross sections, the efficiency threshold is applied to each sample 
separately. Furthermore, since our main motivation for considering scenarios with the second stop 
and the sbottom is to show that the single stop scenario is not overly conservative, we have liberally 
relaxed emin by the square of the branching ratio of Z or W to leptons, for the ^2 and bi samples, 
respectively, in order to avoid artificially penalizing these samples for the smallness of the branching 
ratios into their most distinctive final states. 
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